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Abstract 

The motion of pedestrians is subject to a wide range of influences and exhibits a 
rich phenomenology. To enable precise measurement of the density and velocity 
we use an alternative deflnition using Voronoi diagrams which exhibits smaller 
fluctuations than the standard deflnitions. This method permits examination 
on scales smaller than the pedestrians. We use this method to investigate the 
spatial and temporal variation of the observables at bottlenecks. Experiments 
were performed with 180 test subjects and a wide range of bottleneck parame- 
ters. The anomalous flow through short bottlenecks and non-stationary states 
present with narrow bottlenecks are analysed. 



1. Introduction 

The relation between the flow and the bottleneck width and length and the 
fundamental diagram can be used to validate models and to design facilities. 
Many experiments have been performed to study different aspects of pedestrian 
flow through a bottleneck (for a summary see [1,). Up to now these studies 
have considered bottlenecks with different shape, primarily to resolve ques- 
tions about lane formation, motivation and its effects on the pedestrian flow 
rate EHHIlElElHElElIIOlin]. Statements relating the flow to bottleneck 
parameters should be generally valid and not depend on the initial conditions. 
If a system retains information about its initial conditions then it may exhibit 
a non-stationary state, where quantities such as the density will show a time 
dependence. Analysis of stationarity in these experiments is limited by low par- 
ticipation rates. Effects on stationarity due to the bottleneck width cannot be 
addressed due to the narrow range of bottleneck widths studied. Comparison of 



* Corresponding author 
Email addresses: a. seyfriedOfz- juelich.de (Armin Seyfried), 
b.steffen@fz-juelich.de (Bernhard Steffen), klingsch@uni-wuppertal.de (Wolfram 
Klingsch), rupprechtOuni-wuppertal . de (Tobias Rupprecht), winkens@uni-wuppertal.de 
(Andreas Winkens), in.boltes@fz-juelich.de (Maik Boltes) 



Preprint submitted to Elsevier 



January 20, 2013 



these experiments is also complicated by the differing initial arrangements of the 
pedestrians. Short bottlenecks (to reproduce motion through a doorway) have 
been studied [H [8] and an improved flow rate is seen in agreement with [12 , 
however the former only studied a single bottleneck length so it cannot held to 
be a conclusive result. 

In [12 we observed that the flow exhibited a time dependence (with the 
flow decreasing with time). We were able to make such statements because 
our experiments had a much higher participation level than other experiments 
which have investigated the flow through a bottleneck. Some works have partly 
addressed the stationarity question, through a fitting of an exponential decay to 
the pedestrian velocities |6j or time gaps between successive pedestrians leaving 
the bottleneck [2]. In the case of velocity fits, participation was not great enough 
in order to remove transient behaviour. Due to the massive fluctuations in the 
time gaps no model can be reliably fitted. 

The participants in previous experiments were mixed across social groups. 
Unsurprisingly the mix is not consistent between these experiments, and this 
poses additional problems when trying to compare the experiments. In [S^ the 
effect of population composition on the flow rate through door ways is consid- 
ered. 

If non-stationarity is to be detected a quantity whose fluctuations permit 
the detection of a trend in the data over the course of the experiment is needed. 
Standard definitions of the density exhibit large fluctuations as typically the 
size of the measurement area is comparable with the size of pedestrians and the 
size of the measurement area must be larger than the pedestrians. This limits 
the spatial resolution and stationary states cannot be unambiguously detected. 
An alternative definition using Voronoi diagrams exhibits smaller fluctuations 
exposing non-stationary in our data. This method also permits examination on 
scales smaller than the pedestrians. 

Our experiment was performed in 2006 in the wardroom of the Bergische 
Kaserne Diisseldorf, with a test group of soldiers. The experimental setup al- 
lowed the influence of the bottleneck width and length to be probed. In one 
experiment the width, 6, was varied (from 90 to 250 cm) at a fixed corridor 
length of 400 cm, another experiment investigated varying the corridor length, 
(6,200,400 cm) at a fixed width of 120 cm. Wider bottlenecks with more 
test persons were studied than in previous attempts. The pedestrians begun 
their motion from a waiting area, at a distance d from the bottleneck, to ensure 
an equal initial density for every run, holding areas were marked on the floor 



(dashed regions). The density in this waiting area was 2.6 P/m . Figure 1(b) 



shows a sketch of the experimental setup used to analyze the flow through bot- 
tlenecks, figure p^(a)] shows a still taken from the experiment. The experiments 
investigating the bottleneck width took place on a single day, in order of de- 
creasing bottleneck width. After every 5 experiments a pause was taken. No 
systematic qualitative trends were noted during the course of the experiments. 
The following day the length variation experiments took place in a single run, 
from longest to shortest. 

Our series of experiments have several advantages. The social homogeneity 
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of the test subjects presents a consistent picture of flow characteristics, and the 
high number of participants wih expose any trends in the data. Only one initial 
arrangement of the pedestrians is used. The study of 3 different lengths (6 cm, 
200 cm and 400 cm) allows us to address the problem of stationarity and flow 
through a short bottleneck consistently. The participation level in our exper- 
iment was higher than any previous experiments, combined with the reduced 
fluctuations inherent to integrated method our experiment has the possibility to 
detect non-stationarities in the density. The wider range of bottleneck widths 
studied allow us to address any effects of bottleneck width on stationarity. 



■ -J 

(b) Experimental setup. 
Figure 1: The experimental setup, used in the Bergische Kaserne. 




(a) Combined stills taken from 
experiment. 



2. Definitions 

The density has an underlying distribution, /(x, t), which depends on po- 
sition and time. The exact form of /(x, t) will depend on the measurement 
method and will exhibit fluctuations in time as the system evolves and in space 
on the scale of the individual pedestrians. From the underlying distribution we 
can obtain an average measurement, 

over a measurement area, A, and over a time interval, AT. Varying the size of 
the measurement area or length of measurement time interval, will affect the 
reliability of the density estimate p. For example if the measurement area is 
reduced we must use a longer time period in order to obtain a reliable estimate 
for p, the converse also holds. 

The standard deflnition of the density in a measurement area A, of area 
containing N pedestrians, 

P. = ^, (2) 

exhibits large fluctuations as pedestrians enter and leave the measurement area, 
when the area is of the same order as the size of the pedestrians. These fluc- 
tuations are typically of the same order as the density measurement itself. By 



3 



averaging over time these fluctuations can be removed at the expense of reduced 
time resolution. The underlying density distribution, is a sum of delta 

functions so it can only be applied on scales larger than the pedestrians, limiting 
the potential spatial resolution. 

The focus of this paper is the integrated density method proposed in [13] 
where a density distribution (step function) is assigned to each pedestrian. This 
density distribution is calculated through the Voronoi diagram. At a given time, 
ti, we have a set of positions for each pedestrian {xi(ti), X2(ti), . . . , XMiU)}- We 
compute the Voronoi diagram for these points obtaining a set of cells, A^, for 
each pedestrian i. These cells can be thought of as the personal space belonging 
to each pedestrian. With these cells a density distribution can be deflned for 
each pedestrian. 

The density distribution is the sum of the individual distributions, 

p{x,t) = J2M^'t). (4) 

The density over an area A is deflned as, 

Mt) = (5) 



The integrated method has greatly reduced spatial fluctuations and does not 
require a large AT, to obtain a reliable estimate. The integrated density also 
permits examination on scales much smaller than individual pedestrians, at 
the expense of reduced temporal resolution. In flgure [2] the density is shown, 
calculated using both methods. The integrated density shows much smaller fluc- 
tuations than the classical deflnition. This method provides several advantages. 
The reduced fluctuations mean an instantaneous estimate of the density is pos- 
sible and the presence of non-stationary states can be unambiguously detected, 
something which is not possible with the standard method, see flgure 1. The 
integrated density can also provide microscopic information about the density, 
which will inform the development of microscopic models. 
The pedestrian velocity is deflned as, 

Xi{t ^ At) - Xi{t - At) 
v^{t) = . (6) 

The time difference has been chosen to be At = 0.2 (corresponding to 5 frames in 
advance/behind, for our videos recorded at 25 frames per second). The standard 
velocity is given by 

vsit) = (7) 

Analogous to the integrated density an integrated velocity can be deflned. In- 
dividual contributions to the velocity are 

M^) = l '^f (8) 

^ ^ 10 otherwise ^ ^ 
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(a) b = 120 cm 
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Time / s 

(b) b = 140 cm 



Figure 2: Density computed using standard and integrated definitions. It can be clearly seen 
that the integrated definition fiuctuates less than the standard definition. The non-stationarity 
of the integrated density can be clearly seen (b = 120 cm). The bottleneck is 400 cm long. The 
measurement area is 100 cm x 100 cm, and is centered 100 cm from the bottleneck entrance. 



using the velocity distribution 



(9) 



we can obtain the 



Vv 



v{x)dx 



(10) 



The magnitude of the hydrodynamic flow, 



IJ, 



P\v\, 



(11) 



can also be calculated using the integrated density and velocity. 

Naively calculated Voronoi diagrams possess several features which are in- 
consistent with the concept of a personal space, in particular the cells on the 
border of the diagram. Pedestrians who make up the convex hull have infinite 
cells, other cells maybe large but finite and many cells cross the boundary of 
the room. Three processes are applied to mitigate these problems. By includ- 
ing virtual pedestrians who sit outside the room geometry, all cells are made 
finite. The cells are limited to maximum radius (50 cm) and cut with the room 
geometry. Figure [3] shows a typical diagram after these regulations have been 
applied. 



3. Application 

We detect non-stationarity in the density for the narrowest bottlenecks (90 
cm, 100 cm, 110 cm, 120 cm), which were performed last in the sequence. 
This is an important result, as it points to the emergence of a real physical 
effect, familiarity with the experimental setup does not automatically lead to the 
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Figure 3: Voronoi diagram. The cells are all finite, due to virtual pedestrian regulation, 
cutting with geometry and radius regulation. 



development of a stationary-state. In contrast with [2 where runs were repeated 
to habituate the participants to the experiment, our test subjects exhibited non- 
stationary behaviour after having already got used to the experiment. 

At this point in time we are not studying the transient behaviour of the 
system and this data is removed from consideration. In order to remove transient 
behaviour only data after the first and before the last 30 people are considered. 
This choice of cuts is chosen from observation of the experiments, after the first 
30 have exited, the bottleneck is constantly supplied with pedestrians. During 
the egress of the last 30 people the bottleneck is not supplied by a large reservoir 
of pedestrians. This choice of cuts is preferred as they can be consistently applied 
across all experiments. In figure[4]the density time series are shown using our cut 
criteria. Note that we are discarding more pedestrians than many experiments 
used in total. 
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(a) 6 = 90 cm 
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(c) b = 140 cm 
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(e) b = 200 cm 
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(b) b = 100 cm 
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(f ) b = 250 cm 



Figure 4: Examples of the integrated density timeseries for different bottleneck widths, 6. The 
location of the cuts is shown. 
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3.1. Density, velocity and flow distributions 

With the integrated density method the measurement area may be smaher 
than the pedestrians, something which cannot be done with standard measure- 
ments. By considering the integrated density, velocity and flow calculated over 
small (10 cm x 10 cm) regions and averaging over all stationary state times, 
maps of the observables over the experimental area can be obtained, figures [5] 
and|6] These maps are presented with contours to present additional informa- 
tion, such as the location of the maximum density and directions of steepest 
decline. These maps provide insights into the dynamics of the motion and the 
sensitivity of the integrated quantities to infiuences such as boundary geometry 
and measurement placement. We see that the integrated density is sensitive 
to the geometry, as the Voronoi cells are regulated by the geometry, and the 
cutoff radius applied to the cells. We only present maps for bottlenecks which 
exhibited a stationary state. 

Immediately observable is the strength of the peak in front of the bottleneck. 
Narrower bottlenecks have a higher density in front of the constriction and the 
peak in the density is sharper. The peak in the density extends further into the 
bottleneck for the wider bottlenecks studied. 

In figure [7| the hydrodynamic fiow is shown. The highest fiow is seen in the 
channel and the high fiow region extends further into the waiting area as the 
bottleneck width increases. 
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(c) b = 200 cm (d) b = 250 cm 



jure 6: Integrated velocity maps for bottlenecks of varying width, / = 400 cm. 




(c) b = 200 cm (d) b = 250 cm 

Figure 7: Specific flow, Jg. for bottlenecks of varying width, / = 400 cm. 
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Density and velocity measurements are sensitive to the size and location of 
the measurement area. This is unavoidable, as we are approximating a local 
parameter, e.g. with a measurement over a finite area. Due to this, no 

measure can be considered more correct than any other. When considering 
non-local approximations the dependence on the geometry and location of the 
measurement area must be considered. From the maps it is clear that the size 
and location of the measurement area will infiuence the measured density. 

3.2. Choice of measurement area 

The integrated density follows the same trends as the standard density, how- 
ever we should investigate whether or not any systematic difference can be seen. 
To do this the measurement area, A, must be placed and dimensioned with 
care. Three different effects should be considered in choosing a measurement 
area: the location of the measurement area, the shape of the measurement area 
and the choice of measurement method. The criteria for choosing the location 
of the measurement is simple, we wish to choose an area where the sensitivity 
to measurement location and size is minimized. Due to the symmetry of the 
experimental setup the observation area is chosen to lie directly in front of the 
constriction, symmetrically about the middle. 

We wish to capture information about all the pedestrians passing through 
the bottleneck. As the bottleneck widens the peak of the density spreads in the 
horizontal direction. This implies that we should widen the observation patch 
with the bottleneck. The measurement area is chosen to be 20 cm narrower 
than the bottleneck, and to be 100 cm deep, approximately the depth of two 
pedestrians. 

We wish to place the centre of the measurement area at or near to the peak 
in the density, it is at this point where the measurement will be least sensitive 
to changes in the measurement location. The density maps show that this peak 
lies about 125 cm from the bottleneck entrance. The peak gets a little closer 
(around 100 cm) for the widest bottlenecks. 

In figure [8] the sensitivity of the density to the location of the measurement 
area is shown. The trend of the peak maximum towards the entrance as the 
bottleck widens can be seen. 

In figure [9] the difference between the two measurement methods can be seen, 
using our rectangular measurement area. The integrated density is consistently 
lower than the standard density (although well within fiuctuations). An expla- 
nation lies in the scale sensitivity of the two methods. The integrated method 
is sensitive on scales larger than the measurement area. A pedestrian passing 
near the measurement area will contribute to the integrated density (contribut- 
ing a relatively low density as it sits away from the maximum) but make no 
contribution to the classical density. 

4. Flow through bottleneck 

In [12] it was observed that the fiow through the shortest bottlenecks was 
noticeably greater than the fiow through the longer bottlenecks. This can be 
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Figure 9: Discrepancy between the integrated and standard densities. The error bars represent 
our estimate of the size of the fluctuations. The low density values seen with the narrowest 
bottlenecks are attributed to the non-stationarity seen in these experiments. 
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Figure 10: N(t), the number of people having passed the entrance of the bottleneck. For 
bottlenecks of varying length and equal width (b = 120 cm). 

To inspect the time dependence in a clearer manner we subtract the observed 
pedestrian count N from Jc x where Jc is the average flow over the entire 



experiment. In figure 11 this is plotted. Straight sections correspond to periods 
of constant fiow, positive gradients. For / = 200 cm, 400 cm the horizontal 
sections between t = 20 s and t = 45 s show us that during this period the fiow 
is constant with a value equal to Jc(0- However the / = 6 cm runs exhibit no 
regions of constant fiow. 

The integrated densities ability to probe microscopic scales can be applied 
here to reveal the underlying physical cause. 



In figure [13| the magnitude of the specific fiow, 

\Js\=p\v\, (12) 

is shown for the three varying length experiments performed (calculated over 
10 cm square regions). As is expected for the long bottlenecks the specific fiow 
is constant inside the bottleneck (figures p^3(b)[ p!3(c)] ). In all the experiments 
a protrusion out into the waiting area can be seen. The peak in the specific 
fiow is highest for the shortest bottleneck. These peaks can be seen clearer in 



figure 16 



We now must ask ourselves why this happens. Is the increased fiow due 
to greater pedestrian density at the bottleneck entrance or due to increased 
velocity? In figure [12] the horizontal position where the pedestrians enter the 
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Figure 11: N(t) — JcXt for, the number of people having passed through the bottleneck minus 
a constant flow contribution. The flat and linear areas correspond to regions of stationary 
flow. 



bottleneck is shown. These histograms show that there is no noticeable increase 
in the effective width for the shortest bottleneck. The maps (figures 14 and 15) 
and accompanying slices (figures 17 and 18) demonstrate that the density at the 



entrance to the short bottleneck is significantly lower and the velocity signifi- 
cantly higher. Pedestrians passing through the short bottleneck can side step 
the barrier, in a single pace. This takes them from infront of the bottleneck 
where their motion is restricted by the walls and other pedestrians, to outside 
the bottleneck where their motion is significantly less restricted and they are 
free to move in any direction they choose. This act of side-stepping is respon- 
sible for the lowered density at the entrance to the short bottlenecks, and the 
freer motion responsible for the increased velocity. 
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Figure 12: Crossing points for the individual pedestrians, showing the effective width of the 
corridor. 
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Figure 15: Maps of the integrated velocity. 




Figure 16: Specific fiow measured over 10 cm patches. Slice taken through bottleneck entrance. 




Figure 17: Integrated density measured over 10 cm patches. Slice taken through bottleneck 
entrance. 
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Figure 18: Integrated velocity fiow measured over 10 cm patches. Slice taken through bottle- 
neck entrance. 
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5. Conclusions 

We have shown that the integrated density possess many advantages over 
the classical definition. These are the greatly improved temporal resolution and 
the ability to obtain information about the system over microscopic scales. 

The reduced fluctuations showed us the non-stationarity in the density for 
the narrowest bottlenecks. We do not attribute this non-stationarity to habit- 
uation to the experimental setup, but instead to the emergence of a physical 
effect. This is an important result that merits further investigation, if repro- 
ducible it implies that in many real applications the density is not constant over 
time. The microscopic analysis reveals facets of the density and velocity distri- 
butions which cannot be obtained using standard methods. The density peak 
moves into the bottleneck as the bottleneck widens. The integrated density is 
in agreement with the standard definition, systematic differences between the 
methods can be understood in terms of the scale sensitivity of the methods. 

The anomalous fiow through the short bottlenecks was considered it has 
been shown that the fiow becomes stationary, and that the fiow for the short 
bottleneck is significantly wider. The profile of the specific fiow studied and the 
increased specific fiow was attributed to the presence of a new stepping mode, 
only permitted for short bottlenecks. 
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